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Abstract. Recent observations suggest that gamma ray 
bursts (GRBs) and their afterglows are produced by 
highly relativistic jets emitted in supernova explosions. 
We have proposed that the result of the event is not just 
a compact object plus the ejecta: within days, a fraction 
of the parent star falls back to produce a thick accretion 
disk. The subsequent accretion generates jets and consti- 
tutes the GRB "engine" , as in the observed ejection of rel- 
ativistic "cannonballs" of plasma by microquasars and ac- 
tive galactic nuclei. Here we investigate the production of a 
GRB as the jetted cannonballs exit the supernova shell re- 
heated by their collision with it, emitting highly forward- 
collimated radiation. Each cannonball corresponds to an 
individual pulse in a GRB. We cannot predict the timing 
sequence of these pulses, but the Cannonball Model fares 
very well in describing the total energy, energy spectrum, 
and time-dependence of the individual pulses. 
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1. Introduction 

Once upon a time, Gamma Ray Bursts (GRBs) consti- 
tuted a sheer mystery, whose unassailability was reflected 
in the scores of extremely different ideas proposed to ex- 
plain them. In spite of giant strides in the recent obser- 
vations — the discovery of GRB afterglows (Costa et al. 
1997; van Paradijs et al. 1997), the discovery of the associ- 
ation of GRBs with supernovae (Galama et al. 1998), and 
the measurements of the reedshifts of their host galaxies 
(Metzger et al. 1997) — the origin of GRBs is still an unre- 
solved enigma. In the recent past, the generally accepted 
view has been that GRBs are generated by synchrotron 
emission from fireballs, or firecones, produced by collapses 
or mergers of compact stars (Paczynski 1986; Goodman et 
al. 1987; Meszaros and Rees 1992) by failed supernovae or 
collapsars (Woosley 1993; Woosley and MacFadyen 1999; 
MacFadyen and Woosley 1999, Woosley 1999) or by hy- 
pernova explosions (Paczynski 1998). But various observa- 
tions suggest that most GRBs are produced in supernova 
events by highly collimated ultrarelativistic jets (Shaviv 



and Dar 1995; Dar 1998; Dar and Plaga 1999; Cen 1999; 
Dar and Dc Riijula 2000a and references therein). 

In a previous paper (Dar and De Riijula 2000a) we in- 
troduced a relativistic-cannonball model in which GRBs 
are produced by "cannonballs" (CBs) of baryonic plasma 
emitted subsequently to a core-collapse supernova (SN) 
explosion, and are observable when they happen to point 
close to our direction. There, we concentrated on GRB af- 
terglows — due to bremsstrahlung and synchrotron emis- 
sion from the CBs after they become transparent to their 
own enclosed radiation — to emphasize how, in the case 
of GRB 980425, it might be possible to observe the CBs' 
"superluminal" motion. In this paper we briefly review 
the CB model and we derive its predictions for the prop- 
erties of the 7-rays in a GRB, generated as the forward- 
collimated and blue-shifted thermal radiation from a suc- 
cession of fast-moving, cooling and expanding CBs - 
previously heated by their collision with the SN shell — 
escapes from the transparent outer regions of the shell. 
We study the 7-ray distributions in time, their energy- 
spectrum and the correlations between these two observ- 
ables, showing that the CB model explains the main ob- 
served features of GRBs. 

2. Jets in astrophysics 

Relativistic jets seem to be emitted by all astrophysical 
systems wherein mass is accreted at a high rate from a disk 
onto a central compact object (for a review, see Mirabel 
and Rodriguez 1999a). Highly relativistic jets have been 
observed in galactic sources, such as the microquasars 
GRS 1915+105 (Mirabel and Rodriguez 1994, 1999a,b; 
Rodriguez and Mirabel 1999) and GRO J165-40 (Tingay 
et al. 1995) where mass is accreted onto a stellar black 
hole, and in many active galactic nuclei hosting a massive 
black hole. These jets are not continuous streams: they 
consist of individual blobs of plasma (plasmoids or can- 
nonballs) , and their firing coincides with a sudden removal 
of the accretion-disk material (Belloni 1997; Mirabel and 
Rodriguez 1999b). Cannonballs in microquasars -and pre- 
sumably also in quasars- are emitted in pairs, moving in 
opposite directions. 
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As they travel, microquasar CBs are observed to ex- 
pand at a speed comparable to, or smaller than, the sound 
speed of a relativistic plasma (c/y/3 in their rest system) 
probably because the energy density of their enclosed ra- 
diation is comparable to that of their matter constituency. 
As they become transparent and cool down, the CBs' lat- 
eral size stabilizes to a roughly constant value, presumably 
constrained by magnetic self-containment and/or by the 
ram pressure of the ambient material. Quasar CBs show 
no measurable expansion as they travel, sometimes for as 
long as a million light years (see, e.g., Bridle 2000; Wilson 
et al. 2000). Galactic and quasar CBs expand explosively 
when finally stopped by the material they traverse. 

3. The cannonball model of GRBs 

The ejection of matter in a supernova (SN) explosion is 
not fully understood. The known mechanisms for impart- 
ing the required kinetic energy to the ejecta arc inefficient: 
the theoretical understanding of core-collapse SN events is 
still unsatisfying. It has been proposed (De Riijula 1987; 
Woosley 1993, Dar and De Riijula 2000a and references 
therein) that the result of a SN event is not just a com- 
pact object plus the ejecta: a fraction of the parent star 
may be ejected, but another fraction of its mass may fall 
back onto the newly born compact object. For vanishing 
angular momentum, the free-fall time of a test-particle 
from a parent stellar radius R* onto an object of mass M c 
is: 



tfall 



R2 



-11/2 



; GM r 



1 day 



R* 


3/2 


1.4 Mq~ 


1/2 


10 12 cm_ 




M c 





(1) 



The free-fall time is shorter if the mass of the falling ma- 
terial is not small relative to that of the compact object. 
The fall-time is longer (except for material falling from 
the polar directions) if the specific angular momentum is 
considerably large, as it is in most stars. The estimate of 
Eq.(|l|) is therefore a rough one. 

It is quite natural to suppose that infalling material 
with non- vanishing angular momentum settles into an or- 
biting disk, or a thick torus if its mass is comparable to 
M c . We assume that, as observed in other cases of signif- 
icant accretion onto a compact object (microquasars and 
active galactic nuclei) in which the infalling material is 
processed in a series of "catastrophic" accretions, jets of 
relativistic CBs of plasma are ejected. We presume their 
composition to be "baryonic", as it is in the jets of SS 
433, from which Ly Q and Fe K Q lines have been detected 
(Margon 1984), although the violence of the relativistic 
jetting-process should in our case break most nuclei into 
their constituents. 

The mechanism producing relativistic jets in accre- 
tion processes and its timing-sequence are not understood 



(for suggested possibilities see, e.g., Blandford and Znajek 
1977; Meszaros and Rees 1997). In our model we assume 
that a series of CBs is ejected, each one giving rise to one 
of the "pulses" of a specific GRB. After a few pulses the 
engine runs out of fuel, and the 7-ray activity ceases. The 
timing sequence of the successive pulses we are unable to 
predict, but, as we shall show, the CB model is quite suc- 
cessful in describing the time-dependence of the 7-ray flux 
within single GRB pulses. 

In brief, the CB model is the following. A sequence 
of oppositely-directed pairs of cannonballs is emitted at 
a time tf a ii of 0(1) day after a SN core-collapse. By this 
time the SN outer shell, traveling at a velocity vs ~ c/10 
(see, e.g., Nakamura et al. 2000) has moved to a distance: 



Rs = 2.6 x 10 14 cm 



tfall 
Id" 



10 v s 



(2) 



We adopt Rs = 2.6 x 10 14 cm as a "reference" value, to 
which our results will be scaled. The reference values of 
various relevant parameters — that serve as bench-marks 
to which to scale our results and imply no strong com- 
mitment to their particular choices — are listed in Table I, 
for quick reference. We denote with a bar the actual value 
of a parameter in the units of its reference value so that 
Rs, for instance, means a given SN-shell radius divided by 
2.6 x 10 14 cm. 

Only if traveling at a small angle 9 relative to the line 
of sight, will a CB be visible. As it hits the SN shell, the 
CB slows down and heats up. Its radiation is obscured 
by the shell up to a distance of order one radiation length 
from the shell's outer surface. As this point is reached, the 
GRB is emitted by a CB that continues to travel, expand 
and cool down, its radiation being boosted and collimated 
by the CB's ultrarelativistic motion. We do not discuss 
in this paper the GRB afterglows (Dar and De Riijula 
2000a), the flash of X-ray lines and the achromatic flare 
in the afterglow as the electrons and protons in the GRB 
recombine (Dar and De Riijula 2000b), nor the flux of 
high energy neutrinos and 7-rays produced by the decays 
of pions made in the CB's collision with the SN shell (Dar 
and De Riijula 2000c). 

There are other events in which a variety of GRBs 
could be produced by mechanisms similar to the ones we 
have discussed: large mass accretion episodes in binaries 
including a compact object, mergers of neutron stars with 
neutron stars or black holes (Paczynski 1986, Goodman 
et al. 1987), transitions of neutron stars to hyperon- or 
quark-stars (Dar 1999; Dar and De Riijula, 2000d), etc. 
In each case, the ejected cannonballs would make GRBs 
by hitting stellar winds or envelopes, circumstellar mass or 
light. We discuss only core-collapse SN explosions, as the 
GRBs they would produce by our mechanism, although 
relatively "standard" , satisfactorily reproduce the general 
properties of the heterogeneous ensemble of GRBs. 
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4. Four "clocks" and three energy scales 



Let 7 = 1/yl - f3 2 = E C b/(M C bc 2 ) be the Lorentz factor 
of a CB, that diminishes with time as the CB hits the SN 
shell and as it subsequently plows through the interstellar 
medium. Four clocks ticking at different paces are relevant 
to a CB's history. Let tsN be the local time in the SN rest 
system, tcB the time in the CB's rest system, tob the 
time measured by a nearby observer viewing the CB at an 
angle 9 away from its direction of motion and t the time 
measured by an earthly observer viewing the CB at the 
same angle, but from a "cosmological" distance (redshift 
z ^ 0). Let x be the distance traveled by the CB in the SN 
rest system. The relations between the above quantities 
are: 



dx 

dt SN = 7dt CB = — ; 

(3 c 

dt C B = 5 dt b ; dt = (1 + z) dt b 
where the Doppler factor S is: 
1 27 



5 = 



7 (1-/3 cos 9) (l + 6» 2 7 2 ) 



(3) 



(4) 



and its approximate expression is valid for 9 <C 1 and 
7> 1, the domain of interest here. Notice that for large 7 
and #7 not large, there is an enormous "relativistic aber- 
ration": dt ~ dtsN/7 2 and the observer sees a long CB 
story as a film in extremely fast motion. 

The energy of the photons radiated by a CB in its rest 
system, E CB , their energy in the direction 9 in the local 
SN system, Eg N , and the photon energy, E, measured by 
a cosmologically distant observer, are related by: 



E CB 



E 7 



E 7 



(l + z)E 



(5) 



with S as in Eq.(|). 



5. The making of a GRB 

5.1. Jet energy and CB mass 

Let "jet" stand for the ensemble of CBs emitted in one 
direction in a SN event. If a momentum imbalance be- 
tween the opposite-direction jets is responsible for the 
large peculiar velocities vns ~ 450 ± 90 km s _1 (Lyne and 
Lorimer 1994) of neutron stars born in SNe, the jet kinetic 
energy Ej et must be, as we shall assume for our GRB en- 
gine, larger than ~ 10 52 erg (e.g. Dar and Plaga 1999). 
The jet-emitting process may be "up-down" symmetric to 
a very good approximation, in which case the jet ener- 
gies may be much bigger. There is evidence that in the 
accretion of matter by black holes in quasars (Celotti at 
al. 1997; Ghisellini 2000) and microquasars (Mirabel and 
Rodriguez 1999a, b) the efficiency for the conversion of 
gravitational binding energy into jet energy is surprisingly 
large. If in the production of CBs the central compact ob- 
ject in a SN ingurgitates several solar masses, it is not out 



of the question that Ej ct be as large as M c 2 ~ 1.8 x 10 54 
erg. We shall adopt here a compromise value, 10 53 ergs, 
as the reference jet energy. 

Let 7i n be the Lorentz factor of a cannonball as it is 
fired. Let Ecb = fEjct be the energy of a CB; on av- 
erage GRBs have some five to ten significant pulses, so 
that the fraction f may typically be 1/5 or 1/10. We 
shall adopt Ecb = 10 52 erg as our reference value. For 
this value, the CB's mass is comparable to an Earth 
mass: Mcb ~ 1-8 M 8 (10 3 /7i n ), for a Lorentz factor of 
7in = O(10 3 ), that we shall find to be "typical". 

5.2. CB deceleration by the SN shell 

Let /3i n c be the expansion velocity of a CB, in its rest 
system, as it travels from the point of emission to the 
point at which it reaches the SN shell. For the reference 
value of /3; n , as reported in Table I, we use 1/(10 V^): one 
tenth of the sound velocity of a relativistic plasma. The 
CB reaches the shell with a radius 



Kcb ~ 

7in 



(6) 



^R 2 B X S 



and sweeps up a "target" mass Mt 
M s /3 2 n /(47 2 n ), or some ~ 2.8 x 10" 3 M®, for our refer- 
ence parameter values and ji n = 10 3 . The CB and the 
SN shell are "thick" in the sense of extending over many 
radiation lengths and many nucleon-nucleon interaction 
lengths. A high-energy nucleon suffering successive inter- 
actions in a dilute gas or plasma loses roughly 2/3 of its 
energy to production, with most of the pion energy be- 
ing carried away by the neutrinos in n — ► ji v decays and 
the subsequent /i decays. The electrons from ^ decay and 
the photons from 7r° decay locally deposit roughly 1/3 of 
the original nucleon energy. 

The Lorentz factor of the CB after it has swept the SN 
shell, 7out, is simply the ratio of the total energy to the 
invariant mass (y/s = Mc 2 ) of the outgoing object: 



7out 



Ecb/3 



Ecb/3 



^2M t c 2 E C b/3 + M cb c 4 



(7) 



where we have used Ecb 3> M-pc 2 . Substituting for Mt 
and Mcb as functions of 7i n and (3- ln , one obtains: 



7out ^ 7ii 



2 Ecb 



3/3 2 n M s c 2 + 18E CB 



(8) 



whose limiting values are: 

Tin rt^c-n <1 



7out 



7out 



3 

7b 



(for 6 E CB » At M s c") , 

j_ 

Ecb 



M s 



(for6E CB </3i n M s c 2 ) . (9) 



For our reference 7 out ~ 10 3 , the values of 7; n implied by 
Eqs.(|h may look surprisingly large. But Eqs.(Q) do not 
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depend on Rs : any relativistic jet exiting from the core of a 
SN encounters the same amount of non-collapsed material, 
and must have a 7i n considerably larger than 7 out . 

The very large value of 7i n (~ 3 x 10 3 for our ref- 
erence parameters) implies that the fractional solid an- 
gle covered by a CB as it hits the SN shell is tiny: 
An/(4:7in) ~ lCT 10 , again for our reference parameters. 
This presumably makes it unlikely for successive CBs to 
hit precisely the same spot in the SN shell: CB-CB colli- 
sions and mergers may be the exception, rather than the 
rule. 

5.3. Attenuation of the 7 rays 

The density profile of the outer layers of a SN shell as a 
function of the distance x to the SN center can be mea- 
sured from the photometry, spectroscopy and evolution of 
the SN emissions (see e.g. Nakamura et al. 2000 and ref- 
erences therein). The observations can be fit by a power 
law, x _n , with n ~ 4 to 8. Our results are sensitive to this 
density profile only in the outer region where the SN shell 
becomes transparent (and the measurements are made), 
so that we can adopt the same profile at all x > Rs : 



p(x) = p(Rs)6(x-R s ) 



Rs 



x 



(10) 



The SN-shell grammage still in front of a CB located at x 
is: 

n-l 



Xs(x) 



p(y) dy 



Ms 



4ttR§ 



Rs 



x 



(11) 



For photons in the MeV domain the attenuation length 
is similar, within a factor 2, in all elements from H to Fe 
(Groom et al., 2000), and can be roughly approximated 
by: 



X 7 (E) - 1.0 (E/keV)° M gem 



The value of X 7 (E) in the E = 10 keV to 1 MeV domain 
(2.1 to 9.8 gr/cm 2 ) is close to the attenuation length in a 
hydrogenic plasma (X^ 11 ~ m p /a T ~ 2.6 gr/cm 2 , with m p 
the proton's mass and a T ~ 0.65 x 10 -24 cm 2 the Thomson 
cross-section). Therefore, it makes little difference in prac- 
tice whether or not we take into account that the SN-shell 
material reached by the CB may be ionized by its pre- 
viously emitted radiation. Equating Xg(x) = X 7 (E) and 
solving for x, we define a useful quantity: xt p (E), the posi- 
tion at which the SN shell becomes (one-radiation-length) 
transparent: 



x tp (E) = R s 



M s 1 
4ttR 2 X 7 (E) 



cx E 



-0.33/(n-l) 



(13) 



whose energy dependence is extremely weak. Blue-shifted 
to the SN rest-system, as in Eq.(|B|), GRB photons have 
energies in the MeV range. Let x tp = x tp (l MeV). For 
our reference parameters, some representative results are: 



x tp ~ 2.9 Rs for n = 8, x tp ~ 4.5 Rs for n = 6. At x tp , 
p(x) is orders of magnitude smaller than at x ~ Rs , where 
most of the SN-shell's mass is steeply concentrated. This 
will simplify our discussion, for it is a fair approximation 
to have the CB slow down at heat up close to x = Rs , and 
proceed thereafter unperturbed by the SN-shell material, 
except inasmuch as little of its radiation can escape before 
it reaches x = x tp . At that point, the CB has expanded 
from the radius Rcb of Eq. <M) to a radius at transparency: 



r>t P _ R x t p-Rs x tp -R s 

— *^CB H ~ Pout — ~ Pout 



7out 



7out 



(14) 



some 2.9 x 10 11 cm, for our reference parameters. 
The CB itself becomes transparent to the radiation 
it encloses later, when it reaches a radius R^g — 
[3Mcb OT/(4 7rm p )]5. We expect the CB to stop expand- 
ing at a proper quasi-relativistic rate /3 ou t soon after it be- 
comes transparent and its inner radiation pressure drops 
abruptly: the inertial mildly relativistic transverse motion 
of its matter constituents is slowed-down by interstellar 
material and, perhaps, by self-confining magnetic fields. 

5.4. Total energy of a GRB pulse 

A CB expanding as a quasi-relativistic plasma ought to 
reach the SN shell with a shape (in its rest system) very 
close to spherical. The microscopic description of what 
happens as the CB and the material of the SN shell col- 
lide and coalesce is elaborate (Dar and De Riijula, 2000c). 
Much of the available energy is deposited at the CB's front 
surface by nucleons sharing their energy and 7's from 7r°- 
decay depositing theirs. Electrons from fi decay deposit 
their energy much deeper into the CB. The ionized CB's 
material is hot and dense enough for the deposited energy 
to thermalize very fast. As it impinges the SN shell, a CB 
may have a tendency to get flattened, but the velocities, 

are so sim- 



ile) fj ~ i_ l/(27 2 ), corresponding to 7i n and 70 



ilar that no significant flattening occurs between the time 
the CB hits the shell and the time it reaches the point at 
which the shell becomes transparent to the CB's radiation: 
flattening is subdominant relative to the CB's expansion. 
For the subsequent estimates we approximate the CB as 
a spherical body with a uniform internal temperature. 

The proper temperature To acquired by the CB as 
it hits the SN shell is high enough for the CB's internal- 
radiation energy-density to be much larger than the mass- 
energy density of its matter constituents. Consequently, 
To can be estimated by equating the total internal ra- 
diation energy to the invariant mass in the CB-SN shell 
collision: 



1 7T a 



3 (•'! ftMsc-- 



18E CB ) 3/2 



7put 

R 3 S 



, (15) 



where a 

density constant. The result, T ~ 3.4 keV for our refer- 
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ence parameters, is very insensitive to their values, but for 5.5. Energy and time dependence of a j-ray pulse 



the Rg 3 ^ 4 dependence on the SN-shelPs radius. 

Since the SN-shell's material is highly concentrated 
close to x = Rs, as in Eq.([ll|), we can take To to be 
the temperature at that point. A rough estimate of the 
total energy in a GRB pulse can be obtained as follows. 
While on the part of the shell that is not transparent, 
the CB does not lose much energy via surface radiation, 
so that it expands quasi-adiabatically at roughly con- 
stant RcB(t)T(t). At the point at which the shell be- 
comes transparent, the internal-radiation energy in the 
CB is reduced, from the value y/s of Eq.(^J), to E tp ~ 
y/s Rcb/Rcb- Approximately 1/e of this energy is emit- 
ted thereafter, its value in the CB's rest system is: 



rarest 
pulse 



1 E 



CB 



R-S An Tout 



3 e 7in x tp - Rs Pout Tin 

whose limiting values are: 



(16) 



E£L - 4.5 x 10 45 erg 



3R S 



x tp - Rs 
(for6E CB »AnMs c 2 ), 



1 A, 



E 



£3U * 4.5 x 10 45 erg 



3R S 



x tp - Rs 



Tout Pa 



1 



E 2 



CB 



(for6E CB « An Ms c 2 



ToutAnAout M S 

(17) 



An observer at rest, located at a known luminosity 
distance Dl(z) from the CB and viewing it at an angle 
9 from its direction of motion would measure a "total" 
(time- and energy-integrated) fluence per unit area: 



dF ^ 1+z t 3 

4 ^ D 2 %ulseO 



(18) 



where S = S[j,0] is given, here and in what follows, by 
Eq.(|j) with 7 = 7 out . The "spherical" energy deduced 
from this result would be an overestimate of the true 
energy E™^ by the last factor in Eq.([l8|), which, for 
Tout = 10 3 and #7 ~ O(l), can be as large as ~ 10 9 . 
Enhanced by a factor ranging up to this large number, 
the GRB-pulse energies of Eq.([l6|) can easily reproduce 
the observations, as discussed in detail in Section 7. 

Armed with an expression such as Eq. ( |l8| ) one can em- 
bark in the exercise of studying the extent to which GRBs 
are standard candles, by checking whether the observa- 
tions at fixed redshift are statistically compatible with this 
expression for a uniformly distributed cos 9 distribution. 
But the current number of GRBs with measured redshifts 
is only fifteen, and their deduced total energies are affected 
by absorption, by experimental efficiency and threshold 
effects, etc. We do not, in this paper, attempt such an 
analysis, that has been initiated, with encouraging results, 
by Plaga (2000), who uses — to extract redshifts from a 
large collection of GRBs — the "Cepheid-like" relationship 
between variability and luminosity proposed by Fenimore 
and Ramirez- Ruiz (2000). 



A CB, as it reaches the transparent outskirts of a SN shell, 
is expanding and cooling and its radiation is becoming 
visible to the observer. In what follows it is convenient to 
measure the GRB observer's time, t, setting t = at the 
moment of the encounter of the CB and the SN shell. The 
time of (one-radiation-length) transparency is then: 



Hp 



1 + z x tp - Rs 



7out 5 c 
The CB temperature at t 



Hp 



-t P 



rarest 
pulse 



4™ (Rc P b) 3 



(19) 



(20) 



whose limiting values are: 
3R S 



T 



t P 



-t P 



0.1 keV 



To\ t A! 

Amt 

(for6E CB » A 2 n M s c 2 ), 



Rs 



0.1 keV 



3R S 



"Tout 



E, 



CB 



x tp - R s 
(for6EcB< AnM S c 



An A}ut 
2 u „2\ 



(21) 



The time-dependences of the CB's radius, its temper- 
ature, and the distance x of the CB from the SN's center 
are, for t > 0: 



Rcb [t] — Rcb 
T[t] 

x[t] ~ R s 



Rcb 1 

Hp 



R tp — 

^CB . 

Hp 



Rcb 
tp Rcb [t] ' 
S Tout 



1 



ct 



(22) 



Let the number of photons per unit time and energy, 
assumed to be isotropically emitted by the CB in its rest 
system, be: 



dn 



7 _ 



dE 7 dt C B 



F(EJ B ,T) 



(23) 



Using Eqs.(||||) to change variables to E and t (the j-r&y 
energy and time in the observer's frame), we obtain: 



dn 7 
dE dt 



F E 



: ,T[t] 



(24) 



In the approximation in which the CB's emission in its 
rest system is a thermal distribution from its surface, the 
function F is: 



F(Ej B ,T)c^(R CB [t]) 2 



(E c ^ 2 



Exp{Ej B /T} - 1 



C(3) ^P\^CB/ 
where a = ca/4 is the Stefan-Boltzmann constant. 



(25) 
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The observed energy and time dependence of the pho- 
ton intensity (photon number per unit area, N) of a single 
pulse in a GRB at an angle 9 relative to the CB's motion 
is then predicted to be: 



dN 



1 



dE dt 
dn 



47TD1 



z 2 dn 7 



dEdt ' 
2tt<j [R CB [t] E(l 



z)/S} 2 Abs(E,t) 



dEdt C(3) Exp{E(l + z)/((5T[t])} - 1 
with RcB[t] and T[t] as in Eqs.f^), and where 
Xs(x[t]) 



(26) 
(27) 

(28) 



Ab S (E,t)=E Xp Xt(e(1 + z)) _ 

is the attenuation of the flux in the SN shell. 

For n in Eq.(|l^) as large as the observations indicate 
(n ~ 6), the absorption factor Abs(E 7 , t) rises very sharply 
from to 1 around t — t tp , in which case the width of a 
GRB pulse in energy and time is governed by the expo- 
nential in the denominator of Eq. (p7|) . The argument of 
that exponential can be simply rewritten as Et/H, with: 



H = 



Rs 



C7out 

whose limiting values are: 



(29) 



H- 2.5keVs 



li 



H-2.5keVs 



Tout Pont 

(for6E CB > An M s c 2 
1 E 



CB 



ToutAmtAn 1V1 S 



(30) 



(for6E CB «AnMsc 2 ) . 

6. Some simplifications and approximate 
correlations 



To guide the eye, we give a simplified approximate form of 
Eq.(|l3|), which we do not use in our explicit calculations: 



(Et) 



T Exp{-[t tp /t]- 1 } 0[t]. (31) 



dN 

dEdt " Exp{Et/H} 

The total photon intensity and energy flux are, in this 
approximation: 

dN 



^cxew^Expj-iv/tr 1 } 



F E (t) oc e[t] 



Hp 



Expj-Ittp/tr 1 } 



(32) 



(33) 



Let the peak 7-ray energy at a fixed time during a GRB 
pulse be defined as E^(t) = max [E 2 dl 7 /dEdt]. Its value 
is E^(t) ~ 3.92<5T[t]/(l + z), so that, for t near or after 



Hp- 



E^(t) ^ E 7 (t tp ) e[t] f. 



(34) 



The total "isotropic" energy of a GRB pulse — deduced 
from its observed fluence assuming an isotropic emission — 
can be deduced from Eq. (|l8|), to be: 



4 7^^ 

1 + z 



^pulso 



(35) 



If CBs were "standard candles" with fixed mass, en- 
ergy and velocity of expansion, and if all SN shells had 
the same mass, radius and density distribution, all differ- 
ences between GRB pulses would result from their differ- 
ent distances and angles of observation. For such standard 
candles it follows from Eqs.(||||,|35|) that the observed du- 
rations (half widths at half maximum) of the photon in- 
tensity and of the energy flux density (Ati and Atp), their 
peak values (N p and F p ), and the peak energy (E 7 ) in a 
single GRB pulse are roughly correlated to the total "ob- 
served" isotropic energy (Ej so ) as follows: 



At T oc (1 + z) [E iao \- 
At F oc (1 + z) [E iso ] 



-1/3 



-1/3 



N p oc E iso , 



F p oc [E^] 4 / 3 (1 + z)" 1 , 



K oc [E^] 



1/3 



1 



(36) 
(37) 
(38) 
(39) 
(40) 



These approximate correlations can be tested using the 
sample of 15 GRBs with known redshifts. Because of the 
strong dependence of the CB pulses on the Doppler factor 
and their much weaker dependence on the other parame- 
ters, they may be approximately satisfied (see, e.g. Plaga 
2000) in spite of the fact that CBs and SN shells are likely 
to be sufficiently varied not to result in standard candles. 

Within the standard-candle approximation there is 
also a simple correlation between the rate and the flu- 
ence of GRBs. For the region of the universe that is 
close enough to us to be approximately homogeneous and 
Euclidean, Eq.@ implies that F = E^f S 3 / (4ttD 2 ). 
If CBs were stationary, the corresponding rate of GRB 
pulses with fluence larger than a given F would satisfy 
the well known relation: 



N(> F )~h CB -^ 



E 



pulse 



4ttF 



oc 



-3/2 



(41) 



where hcB is the mean production rate of CBs per unit 
volume. For our highly relativistic CBs, whose "isotropic" 
energy is multiplied by the factor S 3 , Eq.(|4l]) is modified 
to: 



3 2i 



47T 



N(>F )~ - — h CB — 
7 7 Z 



pulse 



4vrF 



oc F, 



-3/2 



(42) 



yielding the same rate-to-fluence relation. The same 
power-law scaling is obtained for the relation between the 
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rate and the peak-energy density-flux from CBs. Both re- 
lations should be approximately satisfied by very bright 
(relatively nearby) GRBs. 

For distant GRBs the above relations arc sensitive to 
the cosmological model, to the not-well-determined SN- 
(or star-formation) rate and to the strong selection ef- 
fect favouring observations of distant GRBs with large 7 
and a small viewing angle 9, and a correspondingly large 
Doppler-enhancement S. Because of this, we do not dis- 
cuss here the rate-to-fluence relation for faint GRBs (Yi, 
1994; Plaga 2000). 

7. Predictions of the Cannonball Model 

Some common properties of GRB pulses (for detailed light 
curves see Kippen 2000; Mallozzi 2000) are observed to be: 

— (a) The GRB fluences, integrated in energy and time, 
lie within one or two orders of magnitude above or 
below 10~ 5 erg/cm 2 (see, e.g., Paciesas et al. 1999). 

— (b) Individual pulses are narrower in time, the higher 
the energy interval of their individual photons (see, 
e.g., Fenimore et al. 1995). 

— (c) Individual pulses rise and peak at earlier time, the 
higher the energy interval of their individual photons 
(see, e.g., Norris et al. 1999; Wu and Fenimore 2000) 

— (d) Individual pulses have smaller photon energies, the 
later the time-interval of observation (see, e.g., Preece 
et al. 1998) . 

— (e) The energy spectrum of GRBs, or of their individ- 
ual pulses, if plotted as E 2 dN/dE, rises with energy 
as E a , with a ~ 1, has a broad peak at E ~ 0.1 to 1 
MeV, and decreases thereafter (see, e.g., Preece 2000). 

— (f ) Most GRBs consist of pulses whose time-behaviour 
is a fast rise followed by an approximately exponential 
decay: a "FRED" shape. Some GRBs have non-FRED, 
roughly time-symmetric pulses (see e.g., Fenimore et 
al. 1995 and references therein) The overwhelming ma- 
jority of GRBs are either made of FRED or non-FRED 
pulses: there are no GRBs with mixed pulse-shapes. 

All of the above items are properties of the CB model. 

In Fig.(|l|) we illustrate item (a) by plotting the total 
fluence, estimated with use of Eq.([l8|) and varying one 
parameter at a time. Naturally, the highest sensitivity is 
that to the viewing angle 9, followed by that to Ecb and 
z. The remaining itemized GRB properties all follow from 
Eq.(27); items (b,c,d) are even apparent in the simplified 
Eq.(31) for the time and energy dependence of the 7-ray 
flux. In Fig.(||) we illustrate items (b) and (c) by plotting 
Eq.(p7|) at three fixed 7-ray energies, for all parameters 
fixed at the reference values of Table I. In Fig.(|J) we sim- 
ilarly illustrate item (d) in a plot at three different times, 
multiples of the time of shell transparency. Item (e) is illus- 
trated in Fig. (4|) where we plot E 2 dN/dE, obtained by in- 
tegrating Eq. ( 27J) over all times; the figure also reports the 
sensitivity to various parameters, by modifying them, one 



at a time, relative to the reference parameters. In Fig. (|5J) 
we illustrate item (f) by plotting dN/dt, obtained by in- 
tegrating Eq.(^7|) over all energies above 30 keV. Once 
again, we vary reference parameters as in Fig.(||). Red- 
shift not being a free parameter specific to our model, 
we separately illustrate in Fig.(^) the z-dependence of the 
time-integrated and energy-integrated versions of Eq. ( p7[ ) . 

A look at Figs.(||) and @ reveals that, for the param- 
eter ranges explored therein, all the predicted GRB-pulsc 
shapes are FREDs and are relatively short in time (frac- 
tions of a second). Yet, these are not general predictions 
of Eq.(^). It is, for instance, quite conceivable that the 
ejection of a shell in a SN explosion be due to one or var- 
ious CBs emitted immediately after core implosion: the 
shock wave induced by their passage through the outer 
shells of the star would trigger their ejection. In that case 
the outgoing shell would be quite disrupted in the "polar" 
directions in which later CBs would result in a GRB. It 
is also possible that a GRB be due to the passage of CBs 
through material expelled by a parent-star's wind, as op- 
posed to the SN shell. In both cases, the density profile of 
the matter traversed by a GRB may be very different from 
that described by a large index n ~ 4 to 8 in Eq. (|l^) , in- 
dicated by observations of complete SN shells, not of their 
small polar regions. In Fig.(0) we illustrate these points by 
plotting dN/dt for n = 2, 3, with the rest of the parame- 
ters at their reference values, and we also give an example 
with n = 6 and a very large viewing angle 9 = 20/7 ou t- 
All three of these time-profiles are quite symmetrical non- 
FREDs and have durations in the few-second range (it is 
also possible to generate long-duration FREDs, as we shall 
see below in the specific case of GRB 980425). 

In Fig.(||) we plot a GRB with 6 CBs, shot at random 
times in a 1.5 s interval and with random values of Ecb 
within a factor of three of our reference value. All other 
parameters in this figure, but the SN-shell density-profile 
index n, have their reference values: the only difference 
between Fig.(||a) and Fig.fQb) is that n = 8 in the former, 
n = 4 in the latter. These figures illustrate the obvious fact 
that the correspondence between CBs and observed pulses 
need not be biunivocal: a CB produces a GRB pulse, but 
an observed pulse can be due to a superposition of CB 
subpulses. Notice that this is also a way to obtain pulses 
that are very wide, or do not have FRED- or symmetrical 
shapes. 

8. Brief comparison to some data 

Comparing a GRB theory with specific GRBs is a tricky 
task, for an obvious reason: GRBs being all different, one 
may be tempted to choose GRBs that fit the theory, rather 
than doing the opposite. In this Section we investigate 
three GRBs with measured redshifts. Of this ensemble, 
we use the highest fluence event (Briggs et al 1999) GRB 
990123 (z = 1.6), to analize the energy distribution; we 
use GRB 980425 (Kippen 2000), whose redshift (Galama 
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et al. 1998) is by far the smallest (z = 0.0085) and yet has 
a conventional fluence (Kippen et al. 1998), to study the 
time-dependence of its single pulse; finally, we use GRB 
990712 (z = 0.4315) to study the correlation between the 
7-ray energy- and time-distributions (Mallozzi 2000) and 
to expose the limitations of the CB model in its present 
simple form. 

8.1. Energy dependence 

The predicted energy spectrum of a GRB is obtained by 
integrating Eq.(^) over all times. The resulting flux dis- 
tribution and the same result weighed with E 2 are com- 
pared with the_GRB 990123 data in_Fig.(|). The param- 
eters used are j3 in = 1/4, (3 out = 1, Ms = 1/5, R s = 1/2, 
E CB = 20, 7 out = 1.5, 6 S = 1.46, n = 6 and z = 1.6. The 
value of Ecb may look large, but this is a multiple-pulse 
GRB and the energy distribution is integrated over all 
pulses: Ecb = 20 corresponds only to twice our reference 
value for Ej et . For these parameters the GRB fluence, as 
estimated via Eq. (|l8|) , is the observed 26.5 10 -5 erg cm~ 2 . 
Since the shape of the energy distribution is insensitive to 
the various parameters, as seen in Fig.(Q), it is easy to 
find many parameter ensembles that result in the same 
prediction: the energy distribution by itself is not a good 
observable to constrain the input, but is, on the other 
hand, a solid test of the model. 

The shape of the energy spectrum dN/dE of Fig.(^a) 
can be easily understood. At the lower energies, the ~1/E 
behaviour is the result of integration over thermal spectra 
with temperatures that decrease with time, see Fig.(||). 
The abrupt decrease of dN/dE at the higher energies re- 
flects the input thermal spectrum at the time and temper- 
ature at which the SN shell starts to become transparent. 

The comparison made in Fig.(^) is quite satisfactory, 
particularly if one realises that many of the higher-energy 
data are but upper limits. In making this figure we used 
the thermal distribution of Eq. ( p5f ) , and the fact that at 
the higher energies the theory may undershoot relative to 
the data is to be expected. Indeed, the CB, in its rest 
system, is subject to a flux of high energy nuclei and 
electrons. While the electrons are being thermalized, they 
should contribute a nonthermal high-energy tail of pho- 
tons emitted via the "free-free" process. Such a power-law 
tail in an otherwise approximately-thermal emission is ob- 
served from young supernova remnants (see, e.g., Dyer et 
al. 2000) and clusters of galaxies (e.g., Fusco-Femiano et 
al. 1999; Rephaeli et al., 1999; Fusco-Femiano et al, 2000), 
both of which are systems wherein a dilute plasma at a 
temperature of 0(1 keV) is exposed to a flux of high en- 
ergy cosmic rays. 

8.2. Time dependence 

In Fig. ([l(]) we compare the single-pulse light curve of GRB 
980425 (Kippen 2000) with the CB theory, obtained by 



integrating Eq.(|27|) over energy, in the 50-300 keV domain. 
The parameters used are (3 in = 1/3, /3 out = 1/2, Mg = 1, 
Rs = 2, Ecb = 10 (corresponding to our reference jet 
energy in a single pulse), 7 out = 1/3, 6>s = 60, n = 8 
and z = 0.0085. For these parameters the GRB fluence, 
as estimated via Eq.(|l8|), is the observed 0.44 x 10~ 5 erg 
cm" 2 . Notice that the value used for the viewing angle 8 
is very large: this is the explanation (Dar and De Riijula 
2000a) why this particular GRB has a normal fluence, in 
spite of how close its progenitor (SN 1998bw) is to us. 

The comparison made in Fig.(|l^) is entirely satisfac- 
tory. The parameter domain giving rise to a light curve 
with a particular shape, height and width is much more 
restricted than the corresponding domain for an energy 
distribution. Yet, we cannot entirely trust the approxi- 
mate parameter values thus extracted, for the reasons to 
be discussed in the next two subsections. 



8.3. The time- energy correlation 

In Fig.([TT|) we compare the single-pulse light curves of 
GRB 990712 (Mallozzi 2000) with the CB theory (the 
continuous red curves), obtained by integrating Eq.(p7|) 
over energy, in the same domains as the data: 20-50 keV 
(BATSE channel 1.1), 50-100 keV (2.2), 100-300 keV (3.3) 
and > 300 keV (4.4). The parameters used are f3 in = 1, 
^ out = 1/3, M s = 1/4, R s = 3, E CB = 50 (corresponding 
to five times our reference jet energy in a single pulse), 
7 out = 1/5, S = 1/2, n = 3, and z = 0.4315. 

The comparison made in Fig.(|ll]) is rather unsatisfac- 
tory, in that the correlation between energy-interval and 
pulse-width is weaker in the observations than it is in the 
predictions. The theoretical correlation, for a thermal in- 
put spectrum, is roughly that implied by the simplified 
expression Eq.(31), that is At AE ~ H, in an obvious no- 
tation. The dashed blue curves in Fig. (|TT|) correspond to 
a modified input in which we have assumed that the CB 
cooling (as discussed in Section 8.5) may be not be linear 
in time, but closer to quadratic, so that At oc l/VAE. 
This modification goes in the right direction, but it is still 
not entirely satisfactory. We have not yet investigated in 
detail how a deviation from an input thermal spectrum 
at high energies -that we discussed in Section 8.1 in com- 
menting Fig.(^)- affects the time-energy correlation. But, 
since a non-thermal high-energy tail broadens the energy- 
distribution at all times, it ought to weaken even further 
the time-energy correlation, as required. 

We have studied the time-energy correlation for other 
single-pulse GRBs, such as 981022, 981221 and 990102. 
They all have a weaker energy-interval to pulse-width cor- 
relation than our model predicts, though the problem is 
most acute for GRB 990712, that we have thus chosen 
to expose the limitations of the CB model in its current 
formulation. 
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8.4- Lessons from the comparison with data 

We conclude from our study of the general properties of 
GRBs in Section 7, and from the three comparisons with 
data in Section 8, that we may have deliniated the cor- 
rect overall energetics of the collision of the CB with the 
SN shell, but our treatment of the time evolution of the 
processes of heating and cooling is oversimplified. A pos- 
teriori, there are many obvious reasons why this ought to 
be the case: the front of the CB is no doubt at a higher 
temperature than its bulk, since the CB is many collision- 
lengths long and is dominantly heated at the front: only 
muons and their decay electrons -but not photons from n 
decay- heat the bulk. The process is not a sudden heating 
followed by continuous cooling, as we assumed. We have 
included cooling by expansion, but not by emission from 
the CB's surface. We have assumed a constant expansion 
velocity, and not attempted to compute an actual expan- 
sion history from plasma dynamics. The CB may not have 
a constant density, it may even be a discontinuous ball of 
"shrapnel". Etc. etc. 

8.5. An alternative simplified model 

The six general properties of GRBs discussed in Section 
7 ought to be quite independent of the complex details 
of the CB's collision with the SN shell, since they only 
capitalize on the overall energetics and on the fact that, as 
it reaches the transparent outskirts of the SN shell, the CB 
is cooling by radiation and expansion. We illustrate this 
point by sketching an alternative model of CB heating and 
cooling, a simplified "surface" model that is in some sense 
the extreme opposite to the simplified "volume" model we 
have discussed in detail. To lighten the discussion, in all 
numerical results in this chapter we fix the parameters to 
their reference values of Table I. 

In its rest frame, the front surface of the CB is bom- 
barded by the nuclei of the SN shell, which have an energy 
m p c 2 7 ~ 1 TeV per nucleon, roughly 1/3 of which (from 
7T° — ► 77 decays) is converted into these 7-rays within 
X p « m p /<7i„(pp) s» 50 gem - , where Cin(pp) is the 
nucleon-nucleon inelastic cross section. These high energy 
photons initiate electromagnetic cascades that eventually 
convert their energy to thermal energy within the CB. The 
radiation length of high energy 7's in hydrogenic plasma, 
dominated by e + e _ pair production, is X 7C ~ 63 g cm , 
comparable to X p . The radiation length of thermalized 
photons in a hydrogenic plasma is X 7 n w m p /er T w 2.6 g 
cm -2 . 

Assume that the quasi-thermal emission rate from the 
CB, within X^ n from its surface, is in dynamical equilib- 
rium with the fraction of energy deposited by the CB's 
collision with the SN shell in that outer layer. The tem- 
perature of the CB's front is then roughly given by: 



T(x) 



(n-1) X 7 m p c 3 [7(x)] 2 cr in (pp) 



6tTX tp X 7e cr 2 



4 


X 




. x tp. 



where 7(x) is a function that decreases monotonically from 
Tin to 7out- Remarkably, only the Lorentz factor of the 
CBs, but neither their mass nor their energy, appear in the 
above expression, except for the fact that, for the result 
to be correct, they must be large enough for the CB to 
pierce the SN shell and remain relativistic. 

For n = 8 the value of xt p is w 3 Rs, and, for t close to 
tt p or later: 



T(t) ~ 0.16 keV 



ttp 


2 


[7(t)l 


_ t _ 




10 3 



(44) 



At the time of transparency this is quite comparable to 
the result of combining Eqs. (pi]) and (p2|). However, only 
for n = 4 does the temperature decrease approximately 
as 1/t. For n > 4 it diminishes faster than 1/t and for 
n = 8 it decreases faster than 1/t 2 , the "faster" being 
due, in both cases, to the effect of a decreasing 7(f). For 
an exact 1 /t 2 behaviour the pulse width narrows with time 
as At cx E~ - 5 and, as we have also seen in Section 8.3, 
this goes in the direction of improving the predicted time- 
energy correlation. In fact, Fenimore et al. (1995) found, 
from a large sample of GRB pulses, that At cx E~°- 46 . 

The total radiated energy, in the CB rest frame, is 
roughly the thermal energy deposition within one radia- 
tion length from its front surface. After attenuation in the 
SN shell, it reduces to: 



Tprest 
pulse 



0"in (pp)^[R' p B ] 2 X 7 m p c 2 7(t) 



3X 



(45) 



7e ' 



(43) 



where X 7 is the radiation length in the obscuring shell 
averaged over the black body spectrum. For a typical 7- 
ray peak energy of E p ~ 1 MeV in the SN rest frame, 
X 7 ~ 10 gem" 2 . Consequently, the CB's radius at trans- 
parency is R^ P B = 4 x 10 11 cm and E£^ e ~ 3 x 10 45 erg, for 
7(t) ~ 10 3 . This is consistent with the results in Eqs.(Q), 
implying that the predicted fluences in the surface-heating 
and volume-heating models are quite similar. The fact 
that the characteristic temperatures of the volume-heating 
and the surface-heating models around the time of trans- 
parency are also similar means that their predicted GRB 
individual-photon energies are comparable and both in 
agreement with the GRB observations. 



9. Conclusions 

In a previous paper (Dar and De Rujula 2000a) we have ar- 
gued that the CB model provides a very good description 
of GRB afterglows, including those whose light curve is 
seen to rise before it drops, as is the case for GRB 970508. 
There we also contended that, in the case of GRB 980425, 
the model provides a strong motivation for the search of 
the super luminal motion of the afterglow-emitting CB, rel- 
ative to the associated supernova: SN1998bw. This would 
be a decisive signature for highly relativistic cannonballs, 
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as opposed to conically spreading jets. We plan to dis- 
cuss in future work other important signatures of the 
CB model: high energy neutrinos and photons during the 
GRB, flare up and X-ray lines in its early afterglow. 

In this paper we have demonstrated that the CB model 
explains the fluence and energy spectrum of GRBs, as well 
as the characteristic properties of their light curves. The 
detailed heating, expanding and cooling of the CB — as it 
hits and sweeps up the SN shell — we have treated only 
in a simplified fashion. As a consequence, the model in 
its present form does not provide a completely satisfac- 
tory quantitative description of the time-energy correla- 
tion. We may not have completely untied the perduring 
Gordian knot of the GRB conundrum, but we have argued 
that we have sliced it open. 
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Parameter 


Symbol 


Value 


SN-shell's mass 


Ms 


10 M 


SN-shell's radius 


Rs 


2.6 x 10 14 cm 


SN-shell's density index 


n 


8 


Outgoing Lorentz factor 


7out 


10 3 


CB's viewing angle 


e 


10~ 3 


CB's energy 


Ecb 


10 52 erg 


Initial v T /c of expansion 


An 


1/(10 \/3) 


Final v T /c of expansion 


Pout 


1/73 


Redshift 


Z 


1 



Table I. List of the "reference" values of various param- 
eters. In the text a barred parameter means its actual 
value divided by its reference value, so that, for instance, 
Mg = 1/2 means the actual mass of the SN shell is taken 
to be 5 M Q . 
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Fig. 1. — Fluence as a function of various parameters, 
in units of 10 -5 erg/cm 2 . The parameters 9 through 7 out 
are in units of their reference values of Table I, thus the 
barred notation. The index n and the redshift z are not 
rescaled. All parameters not being varied are fixed at their 
reference values of Table I, but for 6, fixed at a "typical" 
= 3. 




Fig. 2. — GRB-pulse shape as a function of time, at 
various fixed 7-energies. 
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Fig. 4. — Time-integrated E 2 dN/dE distributions, illus- 
trating the sensitivity to one parameter at a time. The 
absolute vertical scale is arbitrary, but the relative scales 
are not. All parameters not mentioned in each subfigure 
are kept at the reference values of Table I (but for 
which is fixed at 1/ V3) ■ Notice that the shape of the en- 
ergy distribution is always the same, irrespective of the 
parameter values. 
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(g) 
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Fig. 5. — Time distribution dN/dt of a single GRB pulse, 
integrated for all E > 30 keV. The absolute vertical scale 
is arbitrary, but the relative scales are not. All parameters 
not mentioned in each subfigure are kept at the reference 
values of Table I (but for /3; n , which is fixed at 1/V3). For 
ease of comparison, in some subfigures, a curve has been 
rescaled, e.g. in (b) the result for an input Ecb/10 has 
been multiplied by 20. 
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A GRB with 6 CBs 
n = 8 
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Fig. 6. Illustration of the sensitivity to redshift, z, of 
the energy-and time-distributions in a GRB pulse. 



n = 4 




Fig. 7. Examples of parameter values that give rise 
to non-FRED pulse shapes. The individual vertical scales 
are chosen for ease on comparison. All parameters not 
mentioned are as in Table I. 




Fig. 8. — "Synthetic" GRB light curves, generated by 
shooting six CBs at random in a 1.5 s time-interval, and 
with random values of Ecb within a factor of three of 
our reference value. The only difference between (a) and 
(b) is that n = 8 in (a), while n = 4 in (b). All other 
parameters in this figure have their reference values. The 
figure illustrates how a CB produces a GRB pulse, but a 
GRB-pulse may not correspond to a single CB. 
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Fig. 9. — Comparison of theory and observation 
for the time-integrated energy distributions dN/dE and 
E 2 dN/dE, in the case of GRB 990123. Notice that many 
experimental points at the higher energies are only upper 
limits. 
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Fig. 10. — Comparison of theory and observation for 
the light curve of GRB 980425, in the 50-300 keV energy 
interval. 
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Fig. 11. — Comparison of theory and observations for 
the light curve of GRB 980712, in various BATSE energy 
intervals. 



